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Abstract 


A method for eigenvalue and eigenvector approximate analysis for the 
case of repeated eigenvalues with distinct first derivatives is presented. The 
approximate analysis method developed involves a reparameterization of the 
multivariable structural eigenvalue problem in terms of a single positive- 
valued parameter. The resulting equations yield first-order approximations 
to changes in the eigenvalues and the eigenvectors associated with the re- 
peated eigenvalue problem. This work also presents a numerical technique 
that facilitates the definition of an eigenvector derivative for the case of 
repeated eigenvalues with repeated eigenvalue derivatives (of all orders). Ex- 
amples are given which demonstrate the application of such equations for 
sensitivity and approximate analysis. Emphasis is placed on the applica- 
tion of sensitivity analysis to large-scale structural and controls-structures 
optimization problems. 


Introduction 

Future spacecraft will have increasingly higher 
performance requirements and enhanced capabilities. 
Spacecraft such as Space Station Freedom and the 
Earth Observing Satellites will be subject to a wide 
range of disturbance sources, but these spacecraft 
will still be required to meet very stringent pointing 
specifications. Multiple scanning precision instru- 
ments and large-angle payload articulation (robotics) 
are examples of disturbance sources that are ex- 
pected on proposed future spacecraft. Spacecraft 
performance requirements in the presence of mul- 
tiple disturbances present significant technical chal- 
lenges for spacecraft designers. Researchers are cur- 
rently investigating spacecraft design techniques that 
focus on integrating several disciplines to form a 
common interdisciplinary design environment. Inter- 
disciplinary design offers many advantages over con- 
ventional design techniques. The primary advan- 
tage of interdisciplinary design is that it provides 
the framework by which trade-offs, as required in the 
design process, can be explicitly addressed. Specifi- 
cally, design trade-offs on subsystem components can 
be done in such a way as to obtain an optimal de- 
sign for the entire spacecraft; this process allows the 
performance requirements to be more easily met. 

The integrated design of controlled structures is 
a specific application of interdisciplinary design tech- 
niques. A typical example of integrated design of 
controlled structures is the simultaneous controls- 
structures optimization of large flexible space struc- 
tures. This problem has been addressed from a va- 
riety of perspectives. Belvin and Park (ref. 1) have 
proposed a structural tailoring procedure to increase 
the system performance and simultaneously decrease 
the control effort. Maghami et al. (refs. 2 to 4) 


have developed controls-structures integrated design 
methodologies to reduce the amount of control en- 
ergy required to meet pointing requirements while 
simultaneously decreasing the overall mass of the 
structure and increasing its disturbance rejection ca- 
pabilities. Although these techniques are different in 
form and objective, they are computationally inten- 
sive. Most of these computations involve repetitive 
structural eigensolutions. 

The objective of this paper is to present methods 
for increasing the computational efficiency of struc- 
tural optimization and controls-structures integrated 
design methods. Specifically, this paper introduces 
an eigenvalue/eigenvector first-order approximation 
technique for the real, symmetric structural eigen- 
problem for two cases of repeated eigenvalues, and it 
will demonstrate the advantages of the methodology 
in the application of eigensensitivity analysis to facil- 
itate the simultaneous controls-structures integrated 
design process. The two cases of repeated eigenvalue 
problems to be considered are (1) repeated eigen- 
values with distinct eigenvalue derivatives and (2) re- 
peated eigenvalues with repeated eigenvalue deriva- 
tives of all orders (herein referred to as “infinitely 
repeated eigenvalues”). 

This paper also addresses some of the compu- 
tational issues related to eigensensitivity analysis 
for repeated eigenvalue problems. Computational 
issues are addressed by introducing an alternate 
eigenvector derivative matrix formulation that has 
an interesting characteristic (namely, the elimina- 
tion of the numerical burden related to repeated 
factorization of design- variable-dependent eigen- 
sensitivity matrix equations). This formulation elim- 
inates the need for a matrix factorization of the 
large-order sensitivity equations for every design 



variable per eigenvalue/eigenvector pair, and in- 
stead it requires only “back-substitution” for the 
design- variable-dependent right-hand-side vectors . 

Nomenclature 


A 

bi 

E 

ly 

h 

K 

Ke 

k 

L 

M 

^budget 

M e 

Mtot 

n 

ndev 

E 

r 0 

rms 

x* 

Y 

Plos 

Af 

P 

‘Pi 


element cross-sectional area, in 2 

ith design variable 

Young’s modulus, lbf/in 2 

bending moment of inertia about 
Y-axis, in 4 

bending moment of inertia about 
Y-axis, in 4 

global stiffness matrix 

element stiffness matrix 

design variable linking constraint 

element length, in. 

global mass matrix 

total mass budget, slug 

element mass matrix 

total mass, slug 

number of degrees of freedom 

number of structural design variables 

inner radius of tube beam element, in. 

outer radius of tube beam element, in. 

root mean square 

ith structural eigenvector 

orthonormal transformation matrix 

line-of-sight performance measure 

ith structural eigenvalue 

element mass density, lbf/in. 

ith differentiable eigenvector 


Derivation of Eigenvalue and 
Eigenvector Sensitivity Equations 


Repeated eigenvalues may be further subdivided into 
the following three cases: 

1. Repeated eigenvalues with distinct eigenvalue 
derivatives 

2. Repeated eigenvalues with some degree of re- 
peated eigenvalue derivatives 

3. Infinitely repeated eigenvalues 

The simplest case is the distinct eigenvalue prob- 
lem in which eigenvalues are distinct in the design 
space of concern, as shown in figure 1. Methods for 
computing the derivatives of eigenvalues and eigen- 
vectors for the distinct eigenvalue problem have been 
well developed (ref. 5), and therefore they will not be 
addressed here. 



Design variable 

Figure 1. Distinct eigenvalue problem. 


Repeated eigenvalues with distinct eigenvalue 
derivatives often occur at the optimum solution when 
the fundamental (smallest) eigenvalue of the struc- 
ture is maximized with respect to some design vari- 
able. In this case, the value of the fundamental eigen- 
value gradually increases, and at a later stage in the 
design optimization process, it approaches the value 
of the second eigenvalue, as shown in figure 2. In 
this figure, 6° denotes the optimum design variable 
solution, and Ai < A 2 (by definition). Note that the 
eigenvalues Ai and A 2 are the same at 6° and that 
they are only directionally differentiable. 


In the problem to be considered, sensitivities of 
real (generalized) eigenvalues A and eigenvectors x, 
which lie in the null space of the symmetric matrix 
K(6) — AM(6), are examined with respect to vari- 
ations in design parameters (b). Eigenvalues of this 
real, symmetric, structural eigenproblem may be sep- 
arated into two categories: distinct and repeated. 


The computation of eigenvector derivatives in 
the repeated eigenvalue situation is difficult be- 
cause the eigenvectors are not unique ( i.e . , any lin- 
ear combination of eigenvectors corresponding to a 
repeated eigenvalue is also an eigenvector). Re- 
cent publications, however, such as Ojalvo (ref. 6), 
Mills-Curran (ref. 7), Dailey (ref. 8), Juang et al. 
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Design variable 

Figure 2. Repeated eigenvalues with distinct first derivatives. 


that repeat indefinitely. This case is referred to 
as infinitely repeated eigenvalues. Here, a pair of 
eigenvalues and all orders of eigenvalue derivatives 
are equal within the design space of concern, as 
shown in figure 4. Dailey (ref. 8 ) briefly addresses 
the issue when 9X\/9b — OX- 2 /9b. In light of ref- 
erence 9, this equality should be correctly stated 
as d n X\/db n — d n \ 2 [9b n , where n — oo, by not- 
ing that eigenvector derivatives corresponding to re- 
peated eigenvalue derivatives are nonunique. Dailey 
gives an interpretation of nonunique eigenvector 
derivatives that occur with infinitely repeated eigen- 
values analogous to the interpretation of nonunique 
eigenvectors that occur with simple repeated eigen- 
values. Dailey states that any linear combination of 
eigenvector derivatives corresponding to an infinitely 
repeated eigenvalue is also an eigenvector deriva- 
tive. This definition provides an arbitrary para- 
meter, which results in a valid eigenvector derivative 
(ref. 8 ). 



Design variable 


Figure 3. Repeated eigenvalues with repeated first 
derivatives. 

(ref. 9), and Bernard (ref. 10) have proposed numer- 
ical procedures to reorient the eigenvectors as a set 
of differentiable eigenvectors. 

An extension of the last type of eigenvalue prob- 
lem in which not only the eigenvalues but also their 
derivatives are repeated is discussed in reference 9. 
However, it is assumed that the eigenvalue derivative 
is distinct at some finite order. A typical example is 
shown in figure 3, where Ai = A 2 , d\\/db — 6 X 2 /9b, 
and 9 2 X\[9b 2 7 ^ 9 2 X 2 [ 9 b 2 at b°. 

The limiting case of reference 9 is the situa- 
tion in which repeated eigenvalues have derivatives 
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Design variable 

Figure 4. Infinitely repeated eigenvalues. 

The rigid-body modes of a structure can be con- 
sidered as the infinitely repeated case in which the 
eigenvalues remain zero regardless of design changes. 
Structural symmetry is another example that typi- 
cally results in infinitely repeated eigenvalues. Flexi- 
ble symmetric appendages (such as antennas and so- 
lar array masts), which are found on typical space 
structures, usually exhibit infinitely repeated eigen- 
values. The derivatives of these eigenvalues will re- 
peat indefinitely if design changes are made in such 
a way as to preserve structural symmetry. 

In the following two sections, detailed deriva- 
tions of eigenvalue/eigenvector sensitivity equations 
for two types of eigenvalues are presented. For the 
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case of repeated eigenvalues with distinct first deriva- 
tives, analytically exact eigenvector derivative solu- 
tion techniques are presented. These techniques are 
different in form from the modal expansion method 
(ref. 11) and Nelson’s method (ref. 5). The modal ex- 
pansion method assumes that the eigenvector deriva- 
tives can be fully represented as a linear combina- 
tion of a subset of the eigenvectors. The method 
presented in this paper does not make this same as- 
sumption. Instead, the eigenvector derivatives are 
exactly represented as linear combinations of all the 
eigenvectors. In Nelson’s method (ref. 6), decisions 
must be made regarding which components of the 
solution vector should be constrained. The method 
presented herein uses the orthogonality properties 
of the eigenvectors to form a set of linear equality 
constraints. 

For the case of infinitely repeated eigenvalues, 
numerical assumptions are introduced which allow 
the definition of eigenvector derivatives. Emphasis 
is placed upon developing a unified formulation for 
the matrix equations for both types of repeated be- 
havior. All derivations in this paper are based upon 
finite-element representations of structural systems. 

Repeated Eigenvalues With Distinct First 

Derivatives 

This section outlines the derivations of eigen- 
value/eigenvector sensitivity equations for repeated 
eigenvalue problems with distinct first, derivatives. 
For a more detailed derivation of eigenvalue/eigen- 
vector sensitivity equations, as well as eigenvalue/ 
eigenvector approximate analysis, see Kenny (ref. 12). 
In the derivation that follows, it. is assumed that 
the original eigenvectors form a nondefective set. ( i.e . , 
the original set. of linearly independent, eigenvectors 
completely spans an n dimensional space). 

In the case of repeated eigenvalues, the origi- 
nal eigenvectors are, in general, nondifferentiable. 
Therefore, the methods derived for simple (distinct) 
eigenvalue problems are no longer valid. This dif- 
ficulty can be better explained by investigating the 
differences between simple and repeated eigenvalue 
problems. The first, and most, fundamental differ- 
ence is that the eigenvectors of repeated eigenvalue 
problems have an additional degree of nonuniqueness 
compared with those associated with simple eigen- 
values. This degree of nonuniqueness is because any 
linear combination of the eigenvectors is also a valid 
eigenvector. The second difference is related to the 
rank deficiency of the matrix (K — AM), where K 
an d M are the stiffness and mass matrices, respec- 
tively. If repeated eigenvalues occur with a multi- 
plicity of ra, then the matrix (K — AM) will be rank 


deficient, by ra, whereas it. is deficient, by 1 for the 
simple eigenvalue problem. 

To start, the derivation, assume that and X 2 
are the pair of eigenvectors that are associated with 
the repeated eigenvalue A. The differentiable eigen- 
vectors (whose existence is established in ref. 9) 
which are associated with A can then be represented 
as a linear combination of xi and X 2 as follows: 

4>i — Xy; (* = 1,2) (1) 

where X is an n X 2 matrix with x^ and X 2 as its 
columns, and y ; is a 2 X 1 vector. The differen- 
tiable eigenvector should satisfy the n dimensional 
eigenvalue equation 


(K - AM)0,- = 0 (2) 

The literature (refs. 6 to 10 and 12 and 13) shows that 
the derivatives of the repeated eigenvalues are the 
solution of the following reduced eigenvalue problem: 

(K-^M) yi = ° (* = 1,2) (3) 


where; 


K = X r 





(4) 


and 

M = X r MX (5) 

The notation dXi/Ob and 6 X 2 / db, as used in equa- 
tion (3), is defined as the unique eigenvalue deriva- 
tives of the repeated eigenvalue pair at design 
point. 6°, and it. will be used as such throughout, this 
paper. 

Note that because xi and X 2 are linearly inde- 
pendent. and lie within a nondefective set., they can 
be “forced” to become mass orthonormal by any 
orthogonalization process (e.g., the Gram-Schmidt 
process). Therefore, M can be set. equal to a. 2 X 2 
identity matrix without, loss of generality. Equa- 
tion (3) provides a means to find the eigenvalue 
derivatives and the orthonormal transformation ma- 
trix Y with yi and y- 2 , the eigenvectors of equa- 
tion (3), as its columns. The matrix Y will in turn 
allow the definition of a set. of differentiable eigen- 
vectors. If the eigenvalue derivatives are distinct, 
equation (3) yields two simple eigensolutions: 
{dXi/db, yi) and (dX^/Ob, y 2 ). The differentiable 
eigenvectors <p \ — Xyi and (f >2 — Xy 2 correspond 
to the eigenvalue derivatives OXi/Ob and 6 X 2 /Ob, 
respectively. 
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To find the eigenvector derivative d<f»i/db for the 
situation Ai = A 2 = A, one can start with the equa- 
tion for the eigenvector derivative as given below; this 
equation can be obtained by simply differentiating 
equation ( 2 ) with respect to the design parameter b 
such that 


(K- AM) 


d<t>i 

db 


(9A, fdK 8M\ , 




Many methods are available in the literature (refs. 6 
to 9) to find the solution of equation ( 6 ); however, a 
different procedure is presented here. 

This procedure employs the fact that if a so- 
lution to equation ( 6 ) exists, then the eigenvector 
derivative can be rewritten as a linear combina- 
tion of (n — 2 ) eigenvectors and the products of two 
undetermined constants with their corresponding 
differentiable eigenvectors such that 

^ = Uj + c tl (f> 1 + c l2 4>2 {;> =1,2) (7) 

where 

n 

,i = X! c ^3 ( 8 ) 

3 

Note that the particular solution u, is mass orthog- 
onal with respect to the complementary solution 
(Cj '101 + Ci 2 <f> 2 ) such that 


u f M(c i i 0 i + c i 2 0 2 ) 

II 

0 

II 

(9) 

which can be restated as 



ufM 0 ! = 0 

(*' = 1 . 2 ) 

( 10 ) 

ufM 4> 2 — 0 

(*' = 1 . 2 ) 

(ID 

because c,q and c ; 2 are arbitrary in equation 

(7)- 


(See ref. 12 for the proof regarding the existence of a 
solution to eq. ( 6 ).) 

With the aid of equation (7) and the rela- 
tionship dXi/db — 4>J{9K./db — A dN\./db)4>i, equa- 
tion ( 6 ) can be written as 


r rr f rt K 9M\ 1 

(K - AM)u f = \4>J ( 7 ^ - j +i\ ( M <M 


fdK 9 M\ 

( )" (* = 1 , 2 ) ( 12 ) 


db 


db 


to the theory of Lagrange multipliers (ref. 13), a 
unique solution of equation ( 12 ), subject to the linear 
constraints equations ( 10 ) and ( 11 ), can be found as 
the solution of the following equation: 


(K - AM)iij + + iA.i 2 M4>2 


dK 

db 


dM\ 


( a, A J° 


(*' = 1.2) (13) 


in which u, . m q, and m 2 are unknowns, and /in 
and m 2 are the Lagrange multipliers. Nevertheless, 
because equation (13) is a function of the differen- 
tiable eigenvectors <p \ and <p 2 (which are different for 
each design variable), it has to be factored and solved 
for every design variable. Therefore, solving equa- 
tion (13) in a design optimization environment may 
become computationally expensive. Consequently, 
deriving an alternate equation for the eigenvector 
derivatives is desirable. 

To eliminate the design variable dependencies 
from equation (13), a particular solution u; is sought 
by redefining u, as a linear combination of v,-, where 
* — 1,2. In other words, u; = Vy where V is an 
n X 2 matrix with v,-’s as its columns. With this 
new definition and the definition of </>;, equations ( 10 ) 
to ( 12 ) can be written as 

yfv r MX yi =0 (i — 1,2) (14) 

yfv r MXy 2 =0 (* = 1,2) (15) 

and 

r T /9K 9M\ 1 

(K - AMJVy,- = \4>J (-7 - A — J <fc](MXy f ) 

Because eigenvectors yi and y 2 are linearly indepen- 
dent, it can be proved that equations (14) to (15) are 
equivalent to 


v f Mxj = 0 (* = 1,2) (17) 

v f Mx 2 = 0 (* = 1,2) (18) 

To solve equation (16), it is sufficient to solve equa- 
tion (19) such that 


r T fd k 9 M\ 1 

(K - AMJvj = \ 4 >f (-7 - A — J 4 >j\ (Mx f ) 


r K •, M )x. 

rjb rjb J 


(i = 1,2) 

U = 1.2) 


( 19 ) 


whose solution u, has to satisfy the linear equality 
constraints of equations (10) and (11). According 
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Finally, with the aid of the theory of Lagrange multi- 
pliers, the solution of equation (19), subjected to the 
constraints of equations (17) and (18), is equivalent 
to the solution of the following equation: 


'K - AM 

Mxi 

Mx 2 ' 


xfM 

0 

0 

< 

. xjM 

0 

0 




(i =1,2) (20) 


It is worthwhile mentioning here that the derivative 
of the orthogonal condition between <fii and (f> 2 , 
= 0, should be satisfied by the eigenvector 
derivatives d</>i/db and d(f> 2 /db. The eigenvector 
derivatives should preserve the following identity: 


d<t>i 

db 


M(f> 2 




T d M , N 

^"0&~ 02 = O (27) 


which results in an equality with c ^2 and C 21 as 


^ T dM^ 

c 12 + c 21 + 4>i 2 — 0 


(28) 


Note that the left-hand side is independent of design 
variables and that the above equation can be solved 
without knowing the differentiable eigenvectors in 
advance. Recall the eigenvector derivative is defined 
as 

= Vy 8 + c tl 4> 1 + c l2 4>2 (i - 1,2) (21) 

To continue the derivation of d<fii/db, equations 
for the constants c 8 q and in equation (21) need to 
be determined. This determination can be partially 
done by using the first two normalization conditions 
such that 


0fM0,- = l (* = 1,2) 


whose design derivatives yield 


Cll 


1 ( , T dM 

2 


<t> 1 


( 22 ) 


(23) 


and 

C22 = -^(^^ 2 ) (24) 

As for constants c \2 and C 21 , equations similar to 
those given in references 7 and 8 can be derived by 
taking the second derivative of equation (2) to obtain 
the following equalities: 


C12 


<t>l 


<9 2 K 

~W 


<9Ai 

~db~ 


9M <9 2 M 
~db~~ X db 2 


4> 1 


Thus, once c ^2 is found by using equation (25), 
C 21 can be simply obtained via equation (28); this 
process is computationally more efficient than using 
equation (26). 

Note that, although eigenvalue/eigenvector deriva- 
tives associated with repeated eigenvalues can be de- 
termined by equation (3) and equations (20) to (28), 
they cannot be used directly to estimate the change 
in the eigensolution. This estimation will be the 
subject of discussion in the next section. Ex- 
amples are presented subsequently to demonstrate 
the application of the sensitivity equations derived 
above. 

Infinitely Repeated Eigenvalues 

For any structure, the rigid-body modes can be 
found separately from the rest of flexible modes; 
therefore, the discussion that follows focuses on the 
infinitely repeated eigenvalues among the flexible 
modes. 

Similar to the case studied in the previous sec- 
tion, the eigenvalue equation for the infinitely re- 
peated eigenvalues fails to specify the eigenvector 
uniquely. Any linear combination of the eigen- 
vectors xi and X 2 for the eigenvalue A will also be 
an eigensolution. However, for the case of infinitely 
repeated eigenvalues, the eigenvalue derivatives also 
repeat. Therefore, the “eigenvectors” of 




and 


C2lH0T 


+ 2Vl 


( 9K dM 
\db ~ X ~db~ 


d 2 K f d \ 2 
db 2 ^ db 

f dK <9M ' 
\db~ X ~db~ / 



9M <9 2 M 
~db~~ X db 2 



d A 2 
db 


4> 2 

<9Ai 
" ~db 


(25) 


K-fh M )y,=° 


(* — 1 , 2 ) 


(j — 1, ndev ) 


(29) 


are themselves nonunique (i.e. , any linear combina- 
tion of yi and y ~2 will also satisfy the above re- 
duced eigenvalue problem, and the transformation 
from Xj to 4h is no longer unique). Note that the 
limit on subscript j,ndev, is defined as the number 
of independent design variables. 
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In this study, the eigenvectors themselves are 
assumed differentiable as follows: 


4>i - xi (30) 

0 2 = *2 (31) 

Recall that the relationship between 4h and x, is 
given by 4h — Xyy Therefore, the transformation 
matrix for the above assumption is given as 


= 1 (36) 

= 0 (37) 

as given by equations (23) and (24). Nevertheless, 
the orthogonal condition indicated by equation (37) 
imposes an extra constraint on the eigenvector 
derivatives such that 






iT dM , 

0 ! ^ 2 = 0 


(38) 


[yi|y2] 


1 o 
o 1 


(32) 


Furthermore, it is assumed that the derivative of one 
of the infinitely repeated eigenvectors is orthogonal 
to the other repeated eigenvectors. In other words, 
the constants Cij(i j) in equation (7) are assumed 
to be zero. Therefore, the derivatives of the infinitely 
repeated eigenvectors can be expressed as 

9</>l , /ooA 

-ar = v i+ c n x i (33) 


Substituting equations (33) and (34) for d</>i /db 
and d(f> 2 /db, the preceding equation simply implies 
that 

= O (39) 

It is obvious that equation (39) is satisfied exactly 
if clM /db is proportional to M. For the Earth Point- 
ing System (EPS) studied in this paper, numerical 
results show that the values of </>J (d'M./db )</>2 are 
indeed negligibly small, and equations (33) and (34) 
are valid for the infinitely repeated eigenvectors 
considered. 


d(f> 2 , . 

- v 2 + c 22 x 2 (34) 

where v,(i = 1,2) can be solved by equation (20), 

and cn and C 22 can be determined by differentiating 
the following normalization conditions: 

= 1 (35) 


In summary, it is proposed that the derivatives 
of infinitely repeated eigenvectors be expressed by 
equations (33) and (34), provided that the condition 
of equation (39) is not violated. Note that the above 
computational procedure is by no means rigorous, 
although it performs well for the problems studied 
in this work. Future research is required in this 
regard. 


Stiffness and Mass Matrix Derivatives 

In the eigenvalue/eigenvector sensitivity equations derived previously, all equations require derivatives of the 
stiffness and mass matrices. The two most popular methods are the analytical method and the semianalytical 
method. The analytical method, sometimes referred to as the direct differentiation method, uses explicit 
functional relationships between the design variables and the defining matrices. 


Direct or analytical method. For the case in which the design variables are nonlinearly related to the 
stiffness and mass matrices, a chain rule for partial differentiation is required to determine matrix design 
derivatives. As an example, assume that K and M for beam elements have the following general functional 
form: 


K = K (E,G,A,I y ,I z ,L) 
M = M (p,A,I y ,I z ,L) 


(40) 


where G is the shear modulus. Further assume that the generic beam elements defined above are of a tubular 
cross section and that E , G, p, and L are not functions of design parameters (i.e. , E, G, p, and L are constants). 
The terms K and M are then of the following form: 


(41) 


K = K (A, I) 
M = M(T,7) 


7 



where 


A = A(r 0 ,r t ) 
I - I(r 0 ,r t ) 


( 42 ) 


By the chain rule of partial differentiation, the derivative of K and M, with respect to a general radius term bj 
(where bj can be either r 0 or r,, which is the outer and inner radius of the tube element, respectively), is 


<9K 

db 3 

dM 

db. 


<9K / dA ' 
dA \a& Jy 

dM ( dA 
dA 1 db. 


dK dl\ 
dl ) 

dM ( dl 
dl 1 db. 


( 43 ) 


The second derivatives are 
d 2 K 


db 2 


d 2 M 
db 2 


d (dK\ 
db 3 V 9 A ) 

d /dM\ 
db 3 V 9 A ) 


dA 

<9K 

{ d 2 A^ 


[ a 

roK\ 

db] + 

ar 

\ db ) ) 

) + 

K 



d T 
db 


j 


dA 

dM 

( d 2 A' 

, 

' a 

fdM\ 

db 3 

+ ~dA 

\ db ] t 

r 

db 3 



dl 

db,j 


5K / d 2 l\ 

JTydb 2 ) 

/ d 2 l' 
f W\db 2 / 


( 44 ) 


Because the stiffness and mass matrices are linearly related to (A, I), the following terms are zero: 


d (d K 


db 3 \ dA 


db; 


<9K 

W 


d fd M\ 
db 3 V dA ) 

d ddM\ 
db 3 V dl ) 

The second derivatives then can be rewritten as 

d 2 K dK(d 2 A\ 


db 2 

d 2 M 
db 2 


dA \ db 2 , 

/ d 2 A' 
~dA\dbf / 


= 0 


= 0 


= 0 


= 0 


<9K / d 2 l\ 

w{^) 

aM / a 2 /' 

'W\db 2 / 


( 45 ) 


( 46 ) 


The area and bending inertia relationships for a tube beam element are 


A — 7r ( ri — rj 


I - 


4 4 

7r I — rj 


( 47 ) 


where r 0 and r, are, respectively, the outer and inner radii of an element. The first and second derivatives of 
the area and bending inertia terms, with respect to the inner and outer radii, are then 
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dA 


dA 



— 2ir r 0 

— - 

-27T T'i 

dr 0 


dr, 


d 2 A 


d 2 A 



= 2 7T 

— 

— 2 7T 

dr 2 


dr 2 


dl 

3 

dl 

3 


= 7rr^ 

— 

— 7r rf 

dr 0 

o 

dr. 


d 2 I 


d 2 I 



— 37rr 2 

— 

— 37rr. 2 

dr 2 

0 

dr 2 

% 


( 48 ) 


These area and bending inertia derivative terms can be substituted into equations (43) and (46) to complete 
the expressions for the derivatives of the stiffness and the mass with respect to an inner and outer radius design 
parameter. 


A reduction of the number of independent design variables can be obtained by developing relationships 
between design variables. This practice of design variable reduction is commonly referred to as design variable 
linking. A simple example of design variable linking would be to consider the inner radius of a tube a function 
of the outer radius (i.e., r, = kr 0 , where k is a constant that defines the ratio of the inner radius to the outer 
radius, and it must satisfy 0 < k < 1). The area and the bending inertia then take the following functional 
forms : 


A - A[r 0 ,r,(r 0 )\ 
I - I[r 0 ,r,(r 0 )] 


(49) 


With the above functional forms, expressions for the full derivatives of the mass and stiffness matrices with 
respect to the outer radius of a tube beam element can be obtained. Applying the chain rule to equation (41), 
one may obtain as first derivatives 

dK _dK dK f dr, \ \ 

dr 0 dr 0 dr, \dr 0 ) 

\ (50) 

dM _ <9M dM f dr, \ 
dr 0 dr 0 dr, \dr 0 ) ) 

and as second derivatives 


d 2 K _ <9 2 K <9 2 K f dr, 

dr 2 dr 2 dridr 0 \dr 0 


d 2 K d 2 K f dr, 
dr, dr 0 + dr 2 \ dr 0 


dr^ 5K / d 2 r^ 

dr 0 dri \ dr 2 


d 2 M _ <9 2 M <9 2 M / dri 

dr 2 dr 2 dridr 0 \dr 0 


d 2 M d 2 M f dr, 
dr,dr 0 + dr 2 \dr 0 


dr, 

dr 0 


dM ( d 2 r, 
dri \ dr 2 


(51) 


Because the partial derivatives of the stiffness and the mass (with respect to the inner and outer radii) are 
functions of only the inner and outer radii, respectively, this implies that the cross derivative terms are zero. 
The above expressions for the second derivatives then reduce to 


d 2 K 

d 2 K 

d 2 K( 

dr,\ 2 

dr 2 

dr 2 

+ dr 2 \ 

dr 0 J 

d 2 M 

d 2 M 

d 2 M 

1 

> 

f dr, \ 2 

dr 2 

dr 2 

dr 2 

\dr 0 ) 


(52) 


Although the analytical method is an exact method, it can sometimes be too complicated to implement 
if the design variables have complex nonlinear relationships to the stiffness and the mass. However, if design 


9 



variables are linearly related to the stiffness and the mass, the. analytical derivatives can then be obtained 
simply by setting the design variable to unity or zero; this setting depends on whether first or higher order 
derivatives are required. 


For example, if 


then 


and 


K e 


EA _ EA -| 
L L 

EA EA 
L L J 


5Ke 

dA 


r E 

Ei 


L 

L 

= K e 

E 

E 


1 L 

L J 



A—l 


a n K e 

dA n 


o 

o 



. 4=0 


(n > 1) 


(53) 


(54) 


(55) 


Semianalytical method. The semianalytical 
method (ref. 14), which is commonly referred to 
as the finite-difference method for determining stiff- 
ness and mass matrix derivatives, is a technique? 
that greatly simplifies the implementation of eigen- 
value/eigenvector sensitivity analysis. The essence 
of this method is a replacement of the exact stiffness 
and mass matrix derivatives (eqs. (50) and (52)) by 
linearly approximated derivatives. Therefore, non- 
linear design variables can be treated with almost 
the same ease as linear design variables. An example 
of stiffness and mass matrix derivatives formed by a 
forward finite difference is as follows: 


3K ^ 

AK 

K f 

db 3 " 

^ A b 3 

1 

II 

dM ^ 

AM 

_ M : 

dbj * 

^ A&, 

= ~b' 


K° 

b d 

M 


o 


(56) 


where the superscript 0 indicates the quantity defined 
at 6° (the current design), and the asterisk denotes 
the quantity defined at the new design. Specifically, 
K* and M* are the stiffness and mass matrices 
that are formed by a positive perturbation of the 
independent design variable and K° and M° are 
the nominal stiffness and mass matrices. 

The major disadvantage of the semianalytical 
method is that it is an approximate technique, and 
hence the eigenvalue/eigenvector derivatives are no 
longer exact with respect to the finite-element model. 
Therefore, when considering implementation of the 
semianalytical method, one must weigh the advan- 
tage of simplicity versus the disadvantage of de- 
creased accuracy. In the applications to follow, the 


analytical method was chosen because of the high 
numerical accuracy desired. 

Approximate Eigenvalue/ Eigenvector 

Analysis Techniques 

The majority of the computational cost associ- 
ated with any analysis or optimization procedure 
requiring iterative eigensolutions is related to the, 
eigensolutions. This situation is true for even rel- 
atively small problems (<100 degrees of freedom). 
The undeniable fact is that eigensolutions are ex- 
tremely time consuming and computationally bur- 
densome. To alleviate some of this burden, de- 
velopment of methods to efficiently approximate: 
eigensolutions becomes desirable,, One such method 
is the Taylor’s series expansion in which the new 
eigenvalues/eigenvectors are linearly approximated 
as 


ndev <),o 

j — i 3 

ndev o 0 

x* ( n» *?(<)+ 

j — i 3 


(57) 

(58) 


where b* and 6° are the new and current design 
variables, respectively, and A bj is the change in de- 
sign variable bj. The maximum allowable design 
variable perturbation is governed by the linear- 
ity of the function with respect to that variable. 
In general, information regarding the linearity of 
eigenvalues/eigenvectors is not available prior to the 
selection of design variable perturbation sizes. Con- 
sequently, it is virtually impossible to determine 
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general guidelines for their selection. A conservative 
approach, therefore, must be taken to ensure the ac- 
curacy and the stability of the approximated eigen- 
solutions. Based upon the numerical experience of 
this study, the term conservative refers to design 
variable perturbations of <5 percent of nominal 
values (i.e. , A 6 < 0.056°). 

Approximate analysis of eigenvalues /eigenvectors 
in the presence of repeated eigenvalues is compli- 
cated by the apparent discontinuous nature of the 
eigensolution. As mentioned earlier, repeated eigen- 
values are only directionally differentiable, and there- 
fore the use of equations (57) and (58) must be 
severely restricted. Restrictions on equations (57) 
and (58) are required partly because the reduced 
eigenvalue problem (eq. (3)) fails to indicate the 
correspondence between the eigenvalue derivatives 
and the distinct eigenvalues after the design has 
changed. Specifically, to implement equation (57), 
one must have a priori knowledge of the transfor- 
mation from repeated eigenvalues to distinct eigen- 
values after the design has changed. Once this 
transformation is known, the eigenvalue derivatives 
from equation (3) can be properly ordered to yield 
A^ < A^ < . . . < A)/, as required after the design has 
changed. A multidesign variable implementation of 
equation (58) may not be possible because equa- 
tion (58) requires the eigenvector derivatives in the 

ndev 

summation ^ (dx^/dbj) Abj to be associated with 

3 - 1 

the same eigenvector. As seen earlier, differen- 
tiable eigenvectors are design variable dependent, 
and therefore equation (58) may not be a valid ap- 
proximation. A different set of equations, as sug- 
gested by references 6 and 14, should be derived 
to account for the directional dependencies of the 
eigenvalues and eigenvectors. 

A method to reparameterize the multivalued 
eigenproblem into one that is in terms of a single 
positive- valued design parameter will be introduced 
to avoid these directional dependencies. This repara- 
meterization will effectively eliminate the direc- 
tional dependencies of the original multivalued eigen- 
problem, thus allowing the eigenvalues/eigenvectors 
to be approximated using conventional series expan- 
sion techniques (i.e., a specific form of eqs. (57) 
and (58)). 

Let a repeated eigenvalue problem formulated at 
the current design 6 ° be represented as 

k(&°)x° = A°m(&°)x 0 (59) 


Next, introducing a design change A 6 , the new design 
variable b* is then given as 

6 * = 6 ° + Ab (60) 

which yields an eigenvalue problem as 

K( 6 *)x* = A*M( 6 *)x* (61) 

Introducing an intermediate design variable 
6 (e) = 6 ° + eA 6 , the above two eigenvalue equations 
can be collectively represented by a single equation 
such that 


K( 6 (e))x(e) = A(e)M( 6 (e))x(e) (62) 

where e, a real parameter, ranges from 0 to 1. Note 
that the current and new eigenvalue equations can 
be realized by substituting the value e by 0 and 1 , 
respectively. Letting e, which is always positive, be 
assigned as the only design variable, the sensitivity 
equation derived in the previous section can be ap- 
plied here to find the eigenvalue/eigenvector deriva- 
tives of equation (62) with respect to e at e = 0 . 
Note that equation (62) entertains a pair of repeated 
eigenvalues at e = 0 . 

The eigenvalue derivatives of equation (62) at 
e — 0 are the eigenvalues of the following reduced 
eigenvalue problem: 

(K- 7 t -M)y t - = 0 (63) 


where K = X t (K' - AM')X, M = I (where I is the 
identity matrix), and 7 , equals dXi/db. 

The detailed representations of K/ and M/ are 
given as 

dK.(b° + eAb] 


K' = 


de 


e=0 




8=1 


(64) 


and similarly, 


ndev r>, , 


8 = 1 


(65) 
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The, eigenvector derivative <j> f -, corresponding to the 
eigenvalue derivative 77 , is now defined as 

4>\ - Vy ( + c t i 0 i + c t2 4>2 (i = 1 , 2 ) ( 66 ) 

where here y ; is the eigenvector of equation (63), 
and v r - are the solutions of the following equation : 


'K - AM 

Mxi 

Mx 2 " 

1 

\ V( 1 

xfM 

0 

0 


1 H1 

1 

X 

bO 

£ 

0 

0 


{ d i 2 > 


Numerical Verification of Approximate 
Analysis and Applications to Structural 
Optimization 

Several examples will be presented to verify the 
approximate analysis derived for repeated eigenvalue 
problems. In addition, how sensitivity informa- 
tion can be applied in a design optimization loop 
to approximate eigensolutions will be demonstrated. 
These examples will address issues related to success- 
ful implementation of design sensitivity analysis in a 
design optimization environment. 


' — (K' - AM')x; 'j 

° 

0 ) 


(1 = 1,2) (67) 


Note that the above equation is similar to equa- 
tion (20), in which dK jdb and 3M./db are sub- 
stituted by K/ and M/ (defined by eqs. (64) and (65)). 
Furthermore, the constants c ( q and q 2 (where 
i— 1 , 2 ) in equation ( 66 ) can be obtained by us- 
ing equations (23) and (24) and any two of the 
three equations (25) or (26) and (28). In these 
equations, dK /Ob and dtdL/db should be replaced 
by K/ and M/, and d 2 K [ db 2 and d 2 M./db 2 should 
be replaced by and M w , which are defined as 


K" = 


d 2 K{b° + eAb 


d 2 n 

ndevndev 


e=0 


= EE 


l—l 7 = 1 J 


nb;ib A, " A, ‘ 


( 68 ) 


and 


M" 


ndevndev 


EE 


d 2 K 
dbi dbj 


A bi A bj 


(69) 


The first-order approximation of the new eigen- 
value, at e — 1 , can then be done by using the 
Taylor’s series expansion with respect to e about 
e — 0 as 


A?«A° + 7i (1=1,2) (70) 

where 77 is less than j 2 . Therefore, as a result, 
A^ is always less than A^. Similarly, the first-order 
approximation of eigenvectors can be obtained by 

<^0° + ^ (1=1,2) (71) 

where is the differentiable eigenvector at e = 0 
determined by y;, which is the eigenvector 
corresponding to 77 in equation (63). 


Stepped Cantilever Beam 

The goals of this section are to present a simple 
finite-element model with well-understood character- 
istics and to use this model to validate the formu- 
lation for eigenvalue/eigenvector approximate anal- 
ysis in the presence of repeated eigenvalues. Note 
that a priori knowledge regarding the characteris- 
tics of repeated eigenvalues is rarely available. The 
lack of knowledge of this behavior is particularly true 
when considering complex models with multiple de- 
sign variables. With more than two design variables, 
the transformation from repeated to distinct eigen- 
values becomes difficult, if not impossible, to visual- 
ize. Fortunately, a priori knowledge is not required 
for the success of the methods developed in this pa- 
per. Essentially, all that is required is knowledge of 
the particular design variable changes. 

The first part of this validation example will be 
used to introduce the characteristics of the particular 
model being considered. This introduction will be 
done by carefully monitoring the design variable 
perturbations and examining the resulting change in 
the eigensolutions. To facilitate this, only one design 
variable at a time is allowed to be varied. 

The model used in this study is a stepped canti- 
lever beam, which is shown in figure 5. The finite- 
element model consists of 21 nodes, each with three 
degrees of freedom (two translational and one rota- 
tional). The model has been divided into 20 tubular 
beam elements. The ratios between the inner and 
outer radii for both sections are fixed in this study. 
The elements have been subdivided into two groups 
of design variables. The outer radius of the first 10 el- 
ements is design variable 1 ( 61 ), and the outer radius 
of the remaining elements is design variable 2 (& 2 ). 

The initial design parameters for the beam are 
given in table I. At the initial (nominal) design, 
b\ — 0.22699, b 2 — 0.36259, and the frequency of 
bending mode 3 is equal to longitudinal mode 1. 
Bending mode 3 and longitudinal mode 1 correspond 
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Table I. Stepped Cantilever Beam Initial Design Parameters 


Beam parameters 

Design variable 1 

Design variable 2 

Length, in 

Area, in 2 

Unit weight, lbf/in 

!,■ o ' ' ' ' 

Young’s modulus, lbf/in 2 . . . 

5 

.07082 

.92999 

1.4255 x 10“ 3 
10.300 x 10 3 

5 

.18071 

.92999 

9.2807 x 10“ 3 
10.300 x 10 3 




Figure 6. Variations of eigenvalue 3 versus design variables 1 
and 2. 


to eigenvalues 3 and 4, respectively. In this study, 
eigenvalues are considered repeated if their absolute 
difference is <0.01. 

Figures 6 and 7 individually display how eigen- 
values 3 and 4 vary with respect to progressive 
changes in design variables 1 and 2. The combined 
eigenvalue variations are shown in figure 8. These 
figures clearly show that the slopes ( i.e . , eigenvalue 
derivatives ) are discontinuous across the “ridges” 
where the eigenvalues are repeated. 



Figure 7. Variations of eigenvalue 4 versus design variables 1 
and 2. 



Figure 8. Eigenvalue variations versus design variables 1 
and 2. 

Tables II and III present eigenvectors 3 and 4 
at their nominal design (column 3) and with 
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Table II. Eigenvector 3 With Positive Perturbations 
in Design Variables 1 and 2 


Node 

dof 

X 

* 

X N 

* 

1 

1 

0 

0 

0 

1 

2 

0 

0 

0 

5 

1 

0 

0 

7.590 

5 

2 

37.190 

37.125 

0 

9 

1 

0 

0 

14.802 

9 

2 

20.007 

20.000 

0 

13 

1 

0 

0 

19.307 

13 

2 

-13.803 

-13.807 

0 

17 

1 

0 

0 

20.879 

17 

2 

-7.861 

-7.873 

0 

21 

1 

0 

0 

21.412 

21 

2 

25.418 

25.439 

0 


Table III. Eigenvector 4 With Positive Perturbations 
in Design Variables 1 and 2 


Node 

dof 

X 

* 

* 

X i 2 

1 

1 

0 

0 

0 

1 

2 

0 

0 

0 

5 

1 

7.596 

7.595 

0 

5 

2 

0 

0 

37.217 

9 

1 

14.814 

14.810 

0 

9 

2 

0 

0 

20.000 

13 

1 

19.323 

19.319 

0 

13 

2 

0 

0 

-13.785 

17 

1 

20.899 

20.898 

0 

17 

2 

0 

0 

-7.842 

21 

1 

21.433 

21.433 

0 

21 

2 

0 

0 

25.371 


independent positive perturbations in design vari- 
ables 1 and 2 (columns 4 and 5, respectively). In 
these tables, the column headings labeled x? are 

u i 

components of eigenvectors 3 and 4, which are eval- 
uated when design variable i is perturbed by a pos- 
itive A 6j. Also, the column heading labeled “dof” 
refers to either an X degree of freedom (axial defor- 
mation) denoted by 1 or a Y degree of freedom (lat- 
eral deformation) denoted by 2 in the corresponding 
eigenvector. These tables implicitly show the trans- 
formation that must be performed to define a set of 
differentiable eigenvectors. 

To better understand this transformation, one can 
examine the eigenvectors after independent positive 
perturbations in design variables 1 and 2. Tables II 
and III show that with a positive A&i (with A &2 = 0 
(column 4)), the corresponding differentiable eigen- 


vectors are = Xyi = X 3 and <f > 4 = Xy 2 = X 4 . 
Thus, for this type of design change, switching modes 
is not observed. As indicated in column 4 of ta- 
bles II and III, the newly perturbed eigenvectors 
(3 and 4) remain as bending and longitudinal eigen- 
vectors that have nonzero lateral and axial displace- 
ment components, respectively. However, when a 
positive A &2 (with A&i = 0 (column 5)) is intro- 
duced, the corresponding differentiable eigenvectors 
are = Xyi = X 4 and <f > 4 = Xy 2 = X 3 . This trans- 
formation states that the longitudinal mode becomes 
the bending mode, and the bending mode becomes 
the longitudinal mode (i.e . , an exact switch has taken 
place). As indicated in column 5 of table II, the 
newly perturbed eigenvector 3 now becomes a lon- 
gitudinal mode that shows no lateral displacement. 
The same situation can be observed in column 5 
of table III, in which the newly perturbed eigen- 
vector 4 is switched to a bending mode (i.e., zero 
axial displacement). 

With knowledge of the transformation from re- 
peated eigenvalues to distinct eigenvalues, the so- 
lution of the reduced eigenvalue problem (eq. (63)) 
can be properly ordered (i.e., the correspondence be- 
tween eigenvalue derivatives at the repeated eigen- 
value design point and the distinct eigenvalues after 
the perturbation is known). For a positive perturba- 
tion in design variable 1 , the correct ordering of the 
solution of equation (63) is given as follows: 

^ = 195191.2 ( yi = (l, 0 ) T ) 

and 

^1 = 378836.8 (y 2 = (0, 1) T ) 

Similarly, for a positive A& 2 , 

|^ = -237167.5 ( yi = (0,l) T ) 

and 

^1 = 149976.9 (y 2 = (l, 0 ) T ) 

Note that these unitary diagonal transformation 
matrices, Y = [yi,y 2 ], are special cases of the more 
general situation that occurs when the eigenvectors 
are expressed as linear combinations of the pure 
bending and longitudinal modes. In the above discus- 
sion, design variable perturbations have been limited 
to single- variable positive perturbations. A similar 
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analysis can be performed using negative perturba- 
tions. The results will be exactly opposite (i.e . , a 
negative perturbation in design variable 1 will cause 
the two modes to switch, whereas a negative pertur- 
bation in design variable 2 will not cause the modes 
to switch). 

The characteristics of this problem can be de- 
scribed by examining the ratio of the design vari- 
ables b 2 to b\. When the eigenvalues are equal, 
the ratio 62/^1 has a particular value C r , and when 
changes in either b 2 or b 1 increase C r , the original 
modes will switch. Alternatively, when changes in 62 
or b 1 decrease C r , the modes will not switch. 

It should be stressed that the numerical validation 
results that follow do not rely on the knowledge 
gained from this perturbation exercise. The above 
perturbations were given only to provide the reader 
with a better understanding of the characteristics of 
this problem. 

Numerical validation of the method presented for 
eigenvalue/eigenvector approximate analysis is ac- 
complished by investigating two types of design vari- 
able perturbations. The first perturbation contains 
both design variables in such a way as to increase C r , 
thereby causing the original bending and longitudi- 
nal modes to switch. The second perturbation also 
permits both design variables to vary simultaneously ; 
however, the perturbation is such that the ratio C r 
decreases. Regardless of the design change, the meth- 
ods developed in this paper can be used to accu- 
rately predict first-order changes in both the eigen- 
values and the eigenvectors associated with repeated 
eigenvalues. 

Results are presented in the following format. 
The design variable perturbation is presented first. 
The tables that contain the differentiable eigen- 
vectors and their first-order approximations and 
exact values at the perturbed design are presented 
second. Next, a first-order approximation of the 
eigenvalues and their exact values at the newly 
perturbed design are presented. 

Perturbation results for eigenvectors 3 and 4 are 
presented in tables IV to VII. In these tables, col- 
umn 2 is the differentiable eigenvector, column 3 is 
the first-order approximation of the change in the 
eigenvectors as given in equation ( 66 ), and columns 4 
and 5 give the exact eigenvector obtained by solv- 
ing equation ( 2 ) and the approximate eigenvector as 
given by equation (71), respectively. 


Table IV. First-Order Approximation 
for Eigenvector 3 (Perturbation 1) 


X dof 
at 

node 

< f>° 

<t>' 

(h* 

t exact 

^ approximate! 

2 

1.9140 

-0.09726 x 10“ 2 

1.9128 

1.9130 

6 

9.4505 

-.4420 x 10“ 2 

9.4452 

9.4461 

10 

16.5160 

-.6138 x 10“ 2 

16.5080 

16.5099 

14 

19.8110 

-1.2712 x 10“ 2 

19.7970 

19.7983 

18 

21.1320 

-1.8087 x 10“ 2 

21.1120 

21.1139 


Table V. First-Order Approximation 
for Eigenvector 4 (Perturbation f ) 


Y dof 
at 

node 

< f>° 

<t>' 

(h* 

V^exact 

^ approximated 

2 

4.8653 

1.6680 x 10“ 2 

4.8789 

4.8820 

6 

41.7220 

.1165 

41.8120 

41.8385 

10 

7.0410 

-2.7260 x 10“ 2 

7.0083 

7.0137 

14 

-16.0760 

2.8072 x 10“ 2 

-16.0380 

-16.0479 

18 

-.8723 

1.1553 x 10“ 2 

-.8602 

-.8607 


Table VI. First-Order Approximation 
for Eigenvector 3 (Perturbation 2) 


Y dof 
at 

node 

< f>° 

<t>' 

(h* 

V^exact 

^ approximated 

2 

4.8653 

-1.6684 x 10“ 2 

4.8517 

4.8486 

6 

41.7220 

-.1165 

41.6320 

41.6055 

10 

7.0410 

2.7260 x 10“ 2 

7.0726 

7.0683 

14 

-16.0760 

-2.8072 x 10“ 2 

-16.1150 

-16.1041 

18 

-.8723 

-1.1553 x 10“ 2 

-.8844 

-.8838 


Table VIE First-Order Approximation 
for Eigenvector 4 (Perturbation 2) 


X dof 
at 

node 

< f>° 

<t>' 

(h* 

V^exact 

^ approximated 

2 

1.9140 

0.9725 x 10“ 2 

1.9151 

1.9150 

6 

9.4505 

.4420 x 10“ 2 

9.4559 

9.4549 

10 

16.5160 

.6137 x 10“ 2 

16.5240 

16.5221 

14 

19.8110 

1.2712 x 10“ 2 

19.8260 

19.8237 

18 

21.1320 

1.8087 x 10“ 2 

21.1520 

21.1501 


15 



For perturbation 1 (& 2/&1 > CV), the arbitrarily 
selected perturbations of — 0.1 percent in design vari- 
able 1 and 0.1 percent in design variable 2 were cho- 
sen to satisfy the above constraint 62/^1 > CV as 

/ A&i \ / -2.2699 X 10“ 4 

~~ \A6 2 / ~~ V 3.625959 x 10“ 4 

The results of perturbation 1 and eigenvalues 3 
and 4, respectively, are as follows: 

A 3 = A 0 + A 3 = 53488.06 + (-171.99) 

A 3 = 53316.07 (approximated) 

A 3 - 53316.28 (exact) 

A 4 = A 0 + A 4 = 53488.06 + 10.07 
A 4 = 53498.13 (approximated) 

A 4 = 53497.89 (exact) 

For perturbation 2 (b^/b 4 <C r ), arbitrarily se- 
lected perturbations of 0.1 percent in design vari- 
able 1 and — 0.1 percent in design variable 2 were cho- 
sen to satisfy the above constraint b^/b 4 < C r such 
that 


/ A&i \ / 2.2699 x 10 -4 

~\A b 2 )~\ -3.625959 x 10“ 4 

The results for eigenvalues 3 and 4, respectively, 
are as follows: 

A 3 = A 0 + A 3 = 53488.06 + (-10.07) 

A 3 = 53477.99 (approximated) 

A 3 = 53477.75 (exact) 



Figure 9. Earth Pointing System (EPS) structure. 


a 25-m truss, two antennas (15 m and 7.5 m), and 
two antenna supports. Note that the EPS is a hypo- 
thetical structure that was used for early controls- 
structures-interaction (CSI) studies. The finite- 
element model consists of 95 nodes, 223 Euler-type 
beam elements, and 3 design variables. Each node 
has 6 degrees of freedom (dof) for a total of 570 dof. 
The beam elements have tubular cross sections and 
have been divided into three groups of design vari- 
ables (one for each of the major components). The 
outer radius of each group of elements is considered 
a design variable. A design variable linking has been 
chosen such that the inner radius remains propor- 
tional to the outer radius, thereby eliminating the 
need to consider the inner radius as an independent 
design variable. Specifically, the design variable link- 
ing is of the form ?*; = kr 0 , where k is a constant 
that satisfies 0 < k < 1. The design variables for this 
study are as follows: 

• Design variable 1: outer radius of truss 

elements 

• Design variable 2: outer radius of antenna 

elements 


A 4 = A 0 + A 4 = 53488.06 + 171.99 
A 4 = 53660.05 (approximated) 

A 4 = 53660.27 (exact) 

The results presented above clearly show that the 
methods developed in this paper can be used to 
accurately predict first-order changes in eigenvalues 
and eigenvectors for repeated eigenvalue problems 
with distinct eigenvalue derivatives. 

Earth Pointing System Optimization 

The Earth Pointing System (EPS), shown in fig- 
ure 9, is composed of three major component groups, 


• Design variable 3: outer radius of antenna 
support elements 

The characteristics of the structure at the ini- 
tial design are such that it exhibits local vibration 
modes associated with the 15-m and 7.5-m antennas. 
This localized response is a result of the main truss 
section being very rigid at the initial design. Sym- 
metric substructures that exhibit local response are 
almost always responsible for the occurrence of in- 
finitely repeated eigenvalues. As an example, the 
first 20 frequencies of the EPS are tabulated in 
table VIIE The first six frequencies corresponding 
to rigid-body modes are clearly infinitely repeated 
eigenvalues. Moreover, the three pairs of flexible 
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Table VIII. Tabulation of Frequencies for EPS 


Mode 

Frequency, Hz 

Description 

1 to 6 

0 

Rigid body 

7 

.2422 

Torsion 15-m antenna 

8 

.4056 

^-rotation 15-m antenna 

9 

.5652 

A'-rotation 15-m antenna 

10 

.6562 

Torsion 7.5-m antenna 

11 

.8882 

15-m antenna mode 

12 

.8882 

15-m antenna mode 

13 

1.4377 

^-rotation 7.5-m antenna 

14 

1.5360 

A'-rotation 7.5-m antenna 

15 

1.7762 

15-m antenna mode 

16 

1.7762 

15-m antenna mode 

17 

3.0258 

15-m antenna mode 

18 

3.0258 

15-m antenna mode 

19 

3.5132 

15-m antenna mode 

20 

3.5312 

15-m antenna mode 


modes (11 and 12, 15 and 16, and 17 and 18) are 
repeated eigenvalues caused by structural symmetry 
of the 15-m antenna. If design changes are made in 
such a way as to maintain the symmetry of the 15-m 
antenna, those repeated eigenvalues will remain re- 
peated. Note that the 7.5-m antenna also exhibits in- 
finitely repeated eigenvalues, but they are not shown 
in table VIII because they occur outside the 20-mode 
range that is presented. 

This optimization study will demonstrate the 
computational efficiency of using eigensensitivity for 
eigenvalue approximate analysis in the presence of 
infinitely repeated eigenvalues. The objective of this 
problem is to “place” the first eigenvalue of the 15-m 
antenna as high in the frequency range as possible 
with constraints on the total mass of the structure 
and on the first flexible mode of the structure. The 
idea is to leave the flexible modes of the structure 
basically unchanged and “push” the local antenna 
modes higher in the frequency spectrum to avoid ad- 
verse interactions with the control system. To accom- 
plish this, eigenvalue tracking must be enforced ( i.e. , 
the repeated eigenvalue associated with the first an- 
tenna mode must be continually identified as design 
variables are varied as a result of the optimization 
process). Tracking of the first antenna mode is ac- 
complished by selecting the first pair of eigenvalues 
which exhibits infinitely repeated characteristics. 

Mathematically, the optimization problem can be; 
stated as follows: 

maximize 

A j (first infinitely repeated mode) 


subject to 

Aftot < 72.0 slugs (total structural mass) 

Ay < 3.0 (rad/sec) 2 (first flexible mode) 

Ay > 1.0 (rad/sec) 2 

The initial design has the following conditions: 
An = 31.144 (rad/sec)" (first infinitely repeated mode) 
Aftot = 60.2 slugs 
A 7 = 2.32 ( rad/sec )" 
r 1 = rn = r 3 = 1.0021 in. 

For a description of the software used to conduct 
this optimization study and the specific issues re- 
lated to implementing approximate eigensolutions 
within an optimization loop, see appendix A. The 
optimization is performed using the Automated 
Design Synthesis (ADS) software package. The 
optimization technique selected was a linear, ex- 
tended interior penalty function method with a 
Broyden, Fletcher, Goldfarb, and Shanno (BFGS) 
variable metric method for the unconstrained sub- 
problem, followed by a polynomial interpolation for 
the one-dimensional search. 

The optimal design has the following conditions; 

A 12 = 61.054 (rad/sec)" (first, infinitely repeatedmode) 

M tot= 71 9 slu S s 

A7 = 2.002 (rad/sec/ 

?'i = 0.3930 in. (lower bound) rn = 1.4030in.; 7-3 = 1.0033in. 

The results of this study are presented in ta- 
ble IX. Comparisons are based upon using varying 
levels of allowable design variable perturbations for 
the Taylor series eigensolution approximations. An 
allowable design variable perturbation is one in which 
the new design variables remain within a prescribed 
domain of assumed linearity. Figure 10 presents the 
optimization process and graphically shows how ap- 
proximate eigensolutions are used. In this figure, the 
shaded circles represent the domains of assumed lin- 
earity or “linear” neighborhoods. Any design point 
that falls within the domain of assumed linearity may 
be approximated; points outside this region must 
be solved for exactly. Column 2 of table IX shows 
that all solution methods converged to the same op- 
timal design. Table IX presents a performance com- 
parison of four different techniques. Of these four, 
the O-percent-allowable design variable perturbation 
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Table IX. Performance Comparison of Various Levels of Design Variable Perturbations 


Perturbation technique, 
percent approximation 

Objective at 
optimal design, 
(rad/sec) 2 

Function 

evaluations 

requested 

Function 

approximations 

CPU time,* 
sec 

Reduction in CPU 
time, percent 

0 

61.054 

444 

0 

24 643.6 


.5 

61.054 

447 

120 

19 395.2 

21 

2.0 

61.054 

448 

170 

16 729.6 

32 

5.0 

61.054 

456 

296 

11545.7 

53 


*Sun SPARC 2. 



Figure 10. Optimization with approximate eigensolutions. 


technique was by far the most CPU (central pro- 
cessor unit) intensive; the 0.5-percent-allowable de- 
sign variable perturbation was second. A O-percent- 
allowable design variable perturbation technique is 
one in which all eigensolutions are calculated exactly 
(i.e. , no approximate analysis is allowed because the 
area of assumed linearity is zero). The reduction of 
CPU time followed a logical progression; that is, as 
the allowable design variable perturbation increased, 
the CPU time decreased. Note that the decreased 
CPU time is a direct result of the increased num- 
ber of eigensolution approximations (column 4) com- 
pared with the total number of requested function 
evaluations (column 3). As the allowable perturba- 
tion increases, the number of “new” designs that fall 
within the small neighborhood around 6° increases, 
and therefore more functions can be approximated. 
With every eigensolution approximation, a costly ex- 
act eigensolution is avoided, which in turn makes the 
entire optimization process more efficient. In this 
study, a 5-percent-allowable design variable pertur- 
bation in the eigensolution approximations produced 
the best results. This technique yielded more than 
a 50-percent reduction in the CPU time compared 
with the baseline 0-percent technique. 


Finally, note that the baseline 0-percent technique 
would have been much more costly if it were not for 
the particular eigensolver employed. A subspace it- 
eration technique that uses the previous eigenvectors 
as the trial vectors for the current iteration has been 
employed. This method of using the previous eigen- 
vectors greatly decreases the convergence time of the 
new eigensolution. 

Applications of Approximate Analysis to 
Integrated Design of Phase-Zero CSI 
Evolutionary Model 

In the previous example, eigenvalue /eigenvector 
sensitivity analysis was successfully applied in ways 
to significantly reduce the computational burden of 
structural optimization of large finite-element mod- 
els. The computational burden is greatly magnified 
when considering the integrated design of controlled 
structures. For the integrated design techniques pre- 
sented in references 2 to 4, control gains are added 
as design variables to the set of structural design 
variables to form one large simultaneous controls- 
structures optimization problem. One of the most 
difficult tasks in integrated design, or any optimiza- 
tion problem, is the proper selection of realistic ob- 
jective and constraint functions. This selection or 
learning phase adds an additional loop to the design 
process. Therefore, in a design environment that is 
as computationally intensive as integrated design, ev- 
ery effort must be taken to make the optimization as 
efficient as possible. 

Integrated Design Formulation 

The Phase-Zero CSI Evolutionary Model (CEM) 
(ref. 15) was a laboratory test bed designed and 
constructed at Langley Research Center for exper- 
imental validation of control design methods and in- 
tegrated design methodology as part of the NASA 
Controls-Structures Interaction Technology Program 
(ref. 16). The Phase-Zero CEM, shown in figure 11, 
consists of a 62-bay central truss (with each bay 
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Section 6 


Thruster 

accelerometer 

locations 


Section 2 


Figure 11. Phase-Zero CEM. 


10 in. long), two vertical towers, and two horizontal 
booms. The structure is suspended using two cables, 
as shown. A laser source is mounted at the top of one 
of the towers, and a reflector with a mirrored surface 
is mounted on the other tower. The laser beam is 
reflected by the mirrored surface onto a detector sur- 
face 660 in. above the reflector. Eight proportional, 
bidirectional, gas thrusters provide the input actu- 
ation, while collocated servoaccelerometers provide 
output measurements. The first 20 frequencies of the 
Phase-Zero CEM are tabulated in table X. An open- 
loop damping ratio of 0.5 percent is assumed for all 
the modes. 


Table X. Frequencies of Phase-Zero CEM 


Mode 

Frequency, Hz 

Mode 

Frequency, Hz 

1 

0.122 

11 

3.972 

2 

.126 

12 

4.127 

3 

.173 

13 

4.172 

4 

.680 

14 

5.788 

5 

.704 

15 

6.456 

6 

1.159 

16 

6.568 

7 

1.572 

17 

6.777 

8 

1.676 

18 

7.806 

9 

2.085 

19 

8.732 

10 

3.867 

20 

9.396 


To perform the integrated design, the structure 
was divided into seven sections (three sections in the 


main bus and one section each for the two horizontal 
booms and the two vertical towers). Three structural 
design variables were used in each section (namely, 
effective cross-sectional area of the longerons, the 
battens, and the diagonals), thus making a total of 
21 structural design variables. 

The static (or constant-gain) dissipative con- 
troller, which is used for feedback control, employs 
collocated and compatible actuators and sensors and 
consists of feedbacks of the measured attitude vec- 
tor y p and the attitude rate vector y r using constant, 
positive-definite gain matrices G,, and G r . This con- 
troller is robust in the presence of parametric uncer- 
tainties, unmodelled dynamics, and certain types of 
actuator and sensor nonlinearities (ref. 17). Here, 
two of the eight available actuators were used to gen- 
erate persistent white-noise disturbances, and the re- 
maining six actuators were used for feedback control. 
The static dissipative controller uses a 6 X 6 diagonal 
rate-gain matrix with no position feedback. ( Because 
this system has no zero eigenvalues, position feedback 
is not necessary for asymptotic stability.) Thus, in 
the integrated design with the static dissipative con- 
troller, the total number of design variables was 27 
(21 structural plus 6 control design variables). 

An integrated controls-structures design was ob- 
tained by minimizing the steady-state average con- 
trol power hm |jY where Tr( ) denotes 

the trace of ( ), in the presence of white-noise in- 
put disturbances with unit intensity (i.e., standard 
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Design variable 


Figure 12. Structural design variables (initial and final for 
0-percent case). 



Baseline 

CSI-DESIGN, 

percent 


P~j Baseline CSI-DESIGN 
(0-percent case) 


□ 


Approximate eigensolutions with 
5-percent-allowable design 
variable perturbation 


Approximate eigensolutions with 
20-percent-allowable design 



Figure 14. Effect of eigensolution approximations on solution 
integrity. 


the accelerometer outputs were used in the: control 
design. Upper and lower bounds were also placed on 
the structural design variables for safety and prac- 
ticality concerns. Lower bound values were placed 
on these variables to satisfy structural integrity re- 
quirements against buckling and stress failures. Con- 
versely, upper bound values were placed on these 
variables to accommodate design and fabrication 
limitations. 

Mathematically, the integrated design problem 
can be stated as follows: 


Figure 13. CPU time performance comparisons of 0-percent, 
5-percent, and 20-percent eigensolution approximation 
techniques. 

deviation intensity equals 1 lbf), at actuators Num- 
ber 1 and Number 2. (These actuators are located 
at the end of the main bus nearest the laser tower.) 
A constraint was placed on the steady-state line-of- 
sight mis position error lim |jY \e ^J/l 0 sJ/£ s )]}, 

the laser detector (above the structure), for rea- 
sonable steady-state pointing performance. Addi- 
tionally, the total mass of the structure Mt 0 t was 
constrained to facilitate a fair comparison with the 
Phase-Zero CEM design. The six remaining actu- 
ators and the velocity signals (required for feedback 
by the dissipative controllers) obtained by processing 


minimize 

J s , 1 Ui.{ rr [ c ( m,T )]} 

subject to 
and 

M tot < ^budget 

The controls-structures integrated design results 
were obtained using the ADS software package. The 
solutions were obtained using a linear, extended in- 
terior penalty function method with a BFGS vari- 
able metric method for the unconstrained sub- 
problem, followed by a polynomial interpolation for 
the one-dimensional search. 
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Results of the controls-structures integrated de- 
sign are presented in figures 12 to 14. A per- 
formance comparison is based upon using varying 
levels of allowable structural design variable pertur- 
bations. The perturbations considered in this study 
are 5-percent and 20-percent variations in the effec- 
tive areas of the structural members. Solution times 
obtained using these perturbations are compared 
with the baseline CSI-DESIGN code in figure 13. All 
CPU times were obtained using a Sun SPARC 2. 
Results show reductions in CPU times of 72 and 
82 percent for the 5-percent- and 20-percent-allowable 
design variable perturbations, respectively. The ef- 
fect of structural eigenvalue /eigenvector approximate 
analysis on the objective function for the integrated 
design problem is shown in figure 14. This figure 
indicates that approximate analysis with allowable 
design variable perturbations as large as 20 percent 
produces an optimal objective function that is within 
7 percent of the baseline solution. 

Concluding Remarks 

Detailed theoretical derivations of eigenvalue/ 
eigenvector sensitivity derivatives for two cases of the 
real symmetric structural eigenvalue problem have 
been developed in this paper. The two cases ad- 
dressed were the repeated eigenvalues with distinct 
first eigenvalue derivatives and the infinitely repeated 
eigenvalues. An emphasis was placed upon devel- 
oping efficient unified solution techniques for both 
cases. 

The idea of a differentiable eigenvector was in- 
troduced for the case in which repeated eigen- 
values with distinct first eigenvalue derivatives ex- 
isted. A set of differentiable eigenvectors was 
obtained via linear transformation of the original 
eigenvectors. The transformation matrix and the 
eigenvalue derivatives were obtained by solving a re- 
duced eigenvalue problem. The additional conditions 
required to uniquely define the eigenvector derivative 
were obtained by taking the second derivative of the 
structural eigenvalue problem. 

A numerical technique has been proposed for 
the case of infinitely repeated eigenvalues. In this 
proposed technique, the assumption is made that 
the eigenvector derivative corresponding to the first 
eigenvector is orthogonal to the second eigenvector 
for a repeated eigenvalue with a multiplicity of two. 
An additional assumption was introduced regarding 
the “rotation” that must be performed to obtain a 
set of differentiable eigenvectors. More specifically, 
it was assumed that the original set of eigenvectors 
are themselves differentiable eigenvectors. The accu- 
racy of the proposed technique was contingent upon 


satisfying an additional binormalization condition of 
the eigenvectors with respect to the derivative of the 
mass matrix. Approximate analysis results for the 
Earth Pointing System (EPS) model show that the 
proposed technique performs well for this type of 
problem. 

An eigenvalue/eigenvector approximate analysis 
method was then derived in which the perturbed 
eigenvalue problem was continuously transformed 
to an eigenvalue problem with a single positive- 
valued design variable e. The first-order Taylor’s 
series expansion was then applied to estimate the 
change in the repeated eigenvalues and their cor- 
responding eigenvectors caused by the changes in 
design variables. In this technique, the multi- 
variable structural eigenvalue problems defined in the 
neighborhood of the repeated eigenvalues are repa- 
rameterized in terms of the single parameter e, and 
they range from 0 to 1. Consequently, the design 
space becomes one dimensional and positive. There- 
fore, the reordering of the eigenvectors was solely 
dependent upon the magnitudes of the eigenvalue 
derivatives. This technique has been validated by 
numerical examples. 

The uncontrolled EPS structure was used as an 
example to demonstrate the computational efficiency 
of eigenvalue/eigenvector approximate analysis as 
applied to a structural design optimization environ- 
ment. The formulation was based upon a frequency 
objective function that was subject to a mass con- 
straint and additional frequency constraints. The 
goal of this formulation was to demonstrate the effi- 
ciency of using a first-order approximation scheme in 
the presence of repeated eigenvalues and to demon- 
strate the versatility of eigensensitivity analysis by 
implementing an eigenvalue tracking technique based 
on eigenvalue derivative information. The results of 
this study show that more than a 50-percent reduc- 
tion in CPU time can be obtained using a first-order 
approximation scheme. 

Finally, eigenvalue/eigenvector approximate anal- 
ysis was applied to the integrated design of a typi- 
cal large flexible space structure. In this applica- 
tion, approximate analysis was used to facilitate the 
controls-structures optimization of the Phase-Zero 
CEM. The results show that the computational ef- 
ficiency of the design process can be dramatically 
increased by employing the eigenvalue/eigenvector 
approximate analysis. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
March 23, 1994 
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Appendix A 

Controlled Structures Integrated Design 
Software and Implementation of 
Approximate Analysis 

This appendix consists of two sections. The first 
section describes the software that was developed 
to conduct research in the integrated design of con- 
trolled structures. The second section deals with 
specific implementation issues of incorporating struc- 
tural eigenvalue and eigenvector approximate analy- 
sis within the Automated Design Synthesis (ADS) 
optimization loop. 

Controlled Structures Simulation 
Software (CS^) 

The Controlled Structures Simulation Software 
(CS 3 ) test bed (ref. 18) will be used as a research 
tool to investigate various applications of eigen- 
sensitivity analysis to controls-structures optimiza- 
tion. The baseline version of CS 3 consists of three 
modules. These modules include a finite-element 
modeling and structural analysis module, a control 
synthesis module, and an optimization module. The 
CS 3 test bed was developed to take advantage of ef- 
ficient in-core intermodule data transfer and sparse 
matrix storage techniques. For a detailed description 
of the baseline CS 3 package, see reference 18. The 
baseline CS 3 package has evolved to include many 
new capabilities and increased functionality. Most 
of the improvements to CS 3 have been directed to- 
ward increasing its finite-element modeling and con- 
trol synthesis capabilities. Along with its new capa- 
bilities, the CS 3 has also been given a new name. 
This new name, CSI-DESIGN, is more consistent 
with its focus of controls-structures integrated de- 
sign. An eigenvalue/vector sensitivity analysis soft- 
ware (ESAS) module within CSI-DESIGN has been 
developed to implement the eigensensitivity equa- 
tions derived in this paper. The ESAS, which follows 
the format of the CSI-DESIGN package, has been de- 
veloped to utilize efficient in-core data transfer and 
sparse matrix storage techniques. The ESAS con- 
sists of three sections: a decision section, a sensitiv- 
ity analysis section, and an approximate eigenvalue/ 
eigenvector analysis section. The purpose of the de- 
cision section is to evaluate the current status of the 
eigenvalues by determining if they are distinct, re- 
peated with distinct first eigenvalue derivatives, or 
infinitely repeated. Eigenvalue tracking in a general- 
ized format is not enforced. The sensitivity analysis 
section calculates derivatives of the stiffness and mass 
matrices and solves the system of equations required 
to determine the eigenvector derivatives. The types 


of eigenvalue problems considered are distinct, re- 
peated with distinct first eigenvalue derivatives, and 
infinitely repeated. Options are provided for either 
eigenvalue or both eigenvalue/eigenvector sensitiv- 
ity analysis. The types of design variables include 
structural sizing (tube members only) and actuator 
masses (fixed locations). The approximate eigen- 
value/eigenvector analysis section uses eigenvalue/ 
eigenvector derivatives in a truncated Taylor se- 
ries expression to provide linear approximations of 
the eigenvalues and eigenvectors. The approxi- 
mate eigenvalue/eigenvector section provides an op- 
timization interface so that eigensensitivity can be 
used to improve the efficiency of both structure and 
controls-structures optimization. 

Implementation of Approximate 
Analysis in ADS Optimization Loop 

Automated Design Synthesis (ADS) is a general- 
purpose numerical optimization software package 
(ref. 19). The ADS is composed of three main levels: 
the strategy, the optimizer, and the one-dimensional 
search. In the strategy level, the constrained problem 
is successively reduced to a series of linearized sub- 
problems that can be constrained or unconstrained. 
The unconstrained subproblems can be formed by 
any of a variety of penalty function techniques that 
range from interior/exterior penalty function meth- 
ods to the augmented Lagrange multiplier method. 
The constrained subproblem, conversely, can be re- 
parameterized in either linear or quadratic form. The 
optimizer and the one-dimensional search are then 
used to find the minimum of the subproblem. The 
optimizer and the one-dimensional search also have 
a wide variety of options. 

Usually, many function evaluations are required 
to solve the subproblem before convergence is ob- 
tained. In addition to function evaluations, gradi- 
ent information is required to determine search di- 
rections. The ADS allows gradients to be either user 
supplied or internally calculated by finite differences. 
However, ADS does not provide an automatic utility 
to handle simultaneously both user-supplied and in- 
ternally calculated gradients. In some situations, this 
limitation may present less than optimal efficiency. 
As an example, consider the mass minimization of 
a structure subject to nodal point displacement con- 
straints. User-supplied (analytic) gradients of the ob- 
jective function, namely derivatives of the mass with 
respect to sizing design variables, can be determined 
very easily, but gradients of the nodal point displace- 
ments require extensive program development. If 
these displacement gradients are not available ana- 
lytically, the ADS finite-difference option must be 
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chosen. In this case, ADS requires all gradients to 
be internally calculated by finite differences, and the 
analytic mass derivatives cannot be employed. 

A simple technique has been developed to take ad- 
vantage of the efficiency of approximate analysis and 
to alleviate the problem of only partial analytic gra- 
dient information. The technique is outlined below 
as follows: 

1. Allow ADS to calculate finite-difference gradi- 
ents 

2. Approximate, via the Taylor series method 
described in the section entitled “Approxi- 
mate Eigenvalue/Eigenvector Analysis Tech- 
niques,” the perturbed function values for 
those functions whose analytic gradients are 
available 

3. Perform exact function evaluations for those 
functions whose analytic gradients arq not 
available 

To elaborate on the above procedure, a finite- 
difference gradient is calculated as 


dF-i Fi (k°) 

9b~ b o^ A6 


(Al) 


equation (A2). Therefore, the finite-difference calcu- 
lated derivative has “recovered” the analytic deriva- 
tive. This technique is outlined in the flowchart given 
in figure Al. 

To use approximate analysis in the one- 
dimensional line search, a way must exist to ensure 
that the perturbations are within allowable limits. 
This procedure can be accomplished by either set- 
ting limits on the allowable step size for each uncon- 
strained subproblem or by allowing the optimizer to 
take its normal step and approximating the function 
only when this step size is within a given tolerance. 
In this research, the second approach has been chosen 
as the criteria for approximating the eigensolutions. 
A flowchart of this technique is given in figure A2. 



where the tilde implies that the expression is approxi- 
mate. As described in item 2 above, the function may 
be approximated as 




A b 


(A2) 


6 ° 


Substituting equation (A2) into equation (Al) gives 


F- ( b°) + (OF /db ■ 

1 % l v j I \ u 1 %! 


A b - Fi b 


A b 


d_F\ 

Ob; 


(A3) 


6 ° 


Figure Al. Implementation of partial analytic gradient 
information. 



win re I In exact expression in equation (Al), F{(b), Figure A2. Implementation of approximate analysis in ADS 
has been substituted by an approximation given by loop. 
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